Spin-phonon coupling effects in transition-metal perovskites: 
a DFT+[/ and hybrid-functional study 
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Spin-phonon coupling effects, as reflected in phonon frequency shifts between ferromagnetic (FM) 
and G-type antiferromagnetic (AFM) configurations in cubic CaMnOs, SrMnOa, BaMnOs, LaCrOs, 
LaFeOs and La2(CrFe)06, are investigated using density-functional methods. The calculations are 
carried out both with a hybrid-functional (HSE) approach and with a DFT-I-C/ approach using a U 
that has been fitted to HSE calculations. The phonon frequency shifts obtained in going from the 
FM to the AFM spin configuration agree well with those computed directly from the more accurate 
HSE approach, but are obtained with much less computational effort. We find that in the AMnOa 
materials class with A=Ca, Sr, and Ba, this frequency shift decreases as the A cation radius increases 
for the r phonons, while it increases for R-point phonons. In LaMOa with M=Cr, Fe, and Cr/Fe, 
the phonon frequencies at F decrease as the spin order changes from AFM to FM for LaCrOa and 
LaFeOa, but they increase for the double perovskite La2(CrFe)06. We discuss these results and the 
prospects for bulk and superlattice forms of these materials to be useful as multiferroics. 

PACS numbers: 75.85.-|-t, 75.47.Lx, 63.20.dk, 77.80.bg 



I. INTRODUCTION 

Multiferroic materials are compounds showing coexis- 
tence of two or more ferroic orders, e.g., ferroelectric- 
ity together with some form of magnetic order such as 
ferromagnetism or antiferromagnetism. Presently such 
materials are attracting enormous attention due to their 
potential for advanced device applications and because 
they offer a rich playground from a fundamental physics 
point of view. Possible applications tend to focus on the 
magnetoelectric (ME) coupling, which may pave the way 
for control of the magnetization by an applied electric 
field in spintronic devices ^i*^ although other applications 
such as four-state memories^ are also of interest. 

Although intrinsic multiferroic materials are highly de- 
sirable, they are generally scarce. One often-proposed 
reason may be that partially filled Sd shells favor mag- 
netism, while the best-known ferroelectric (FE) per- 
ovskites have a Sd" configuration for the transition metal 
(e.g. Ti^+, Nb^+, etCi^). However, it has recently 
been shown^ that some magnetic perovskite oxides dis- 
play incipient or actual FE instabilities, clearly indi- 
cating that there are ways around the usually-invoked 
incompatibility between ferroelectricity and magnetism. 
Some of these compounds, namely CaMnOs, SrMnOa 
and BaMnOa, will be considered in this work. 

Even if several microscopic mechanisms have been 
identified for the occurrence of ferroelectricity in mag- 
netic materials alternative routes are needed in or- 
der to optimize materials for functional device applica- 
tion. This could be done cither by focusing on new 
classes of compounds such as hybrid organic-inorganic 
materials^-^-iS or by modifying already-known materials 
to engineer specific properties. In the latter direction, 
an intriguing possibility is to start with non-polar ma- 



terials and then induce multiple non-polar instabilities; 
under appropriate circumstances this can produce a fer- 
roelectric polarization, as first predicted in Ref. [O based 
on general group theory arguments and analyzed in the 
SrBi2Nb2 09 compound by means of a symmetry analysis 
combined with density-functional theory calculations by 
Perez-Mato et alJ^ Here the ferroelectricity was found to 
arise from the interplay of several degrees of freedom, not 
all of them associated with unstable or nearly-unstable 
modes. In particular, a coupling between polarization 
and two octahedral-rotation modes was invoked to ex- 
plain the behavior,!^ Bousquet et al. have demonstrated 
that ferroelectricity is produced by local rotational modes 
in a SrTiOa/PbTiOs superlattice Although in most 
improper ferroelectrics a single primary order parame- 
ter induces the polarization;^ Benedek and Fennie pro- 
posed that the combination of two lattice rotations, nei- 
ther of which produces ferroelectric properties individ- 
ually, can induce a ME coupling, weak ferromagnetism, 
and ferroelectricity..-- Indeed, we now know that rota- 
tions of the oxygen octahedra, in combinationi^^— and 
even individually ) -^^'^^ can produce ferroelectricity, mod- 
ify the magnetic order, and favor magnetoelectricity. 

Another route to creating new multiferroic materials 
may be to exploit the coupling between polarization, 
strain, and spin degrees of freedom. A strong dependence 
of the lowest-frequency polar phonon frequency on epi- 
taxial strain^^ is common in paraelectric (PE) perovskite 
oxides, and can sometimes be exploited to drive the 
system ferroelectric, a phenomenon known as epitaxial- 
strain-induced ferroelectricity.^^ In a magnetic system 
that also has a strong spin-phonon coupling, i.e., a strong 
dependence of the polar phonon frequencies on spin or- 
der, the magnetic order may be capable of tipping the 
balance between PE and FE states. For example, con- 
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sider a system that has two competing ground states, one 
of which is antiferromagnetic (AFM) and PE while the 
other is ferromagnetic (FM) and FE, and assume that the 
spin-phonon couphng is such that the lowest-frequency 
polar phonon is softer for FM ordering than for AFM 
ordering. In such a case, the epitaxial-strain enhance- 
ment of a polar instability may lead to a lowering of the 
energy of the FM-FE state below that of the AFM-PE 
phase. Spin-phonon mechanisms of this kind have been 
powerfully exploited for the design of novel multiferroics 
in Refs. EtS 

Clearly, it is desirable to have a large spin-phonon cou- 
pling, in terms of a large shift Aw of phonon frequen- 
cies upon change of the magnetic order. Interestingly, 
recently computed spin-phonon couplings in ME com- 
pounds seems to be strikingly large. For example, a 
Auj of about 60 cm~^ has been reported for the dou- 
ble perovskite La2NiMn06.^- Furthermore, Alu values of 
about 200 cm~^ have been computed for a cubic phase 
of SrMnOai^ These values appear to be anomalously 
large when compared to phonon splittings across mag- 
netic phase transitions in other oxides, which are in the 
5-30 cm~^ rangei^i^ 

In the search for materials with large spin-phonon cou- 
pling, first-principles based calculations play a prominent 
role since they can pinpoint promising candidates simply 
by inspecting the dependence of the phonon frequencies 
on the magnetic order. However, a serious bottleneck 
appears. To obtain a reliable and accurate description 
of these effects, it is important to describe the struc- 
tural, electronic and magnetic properties on an equal 
footing. This is especially true for transition-metal ox- 
ides, which are the usual target materials for the spin- 
phonon driven ferroelectric- ferromagnet. In this case, it 
is well known that the localized nature of the 3d elec- 
tronic states, or loosely speaking, the "correlated" nature 
of these compounds, poses serious limits to the appli- 
cability of common density-functional methods like the 
local density approximation (LDA) or generalized gradi- 
ent approximation (GGA). In fact, these standard ap- 
proximations introduce a spurious Coulomb interaction 
of the electron with its own charge, i.e., the electrostatic 
self-interaction is not entirely compensated. This causes 
fairly large errors for localized states (e.g., Mn d states). 
It tends to destabilize the orbitals and decreases their 
binding energy, leading to an overdelocalization of the 
charge density.'^^, 

One common way out is the use of the DFT+J7 
method,^"* where a Hubbard-like U term is introduced 
into the DFT energy functional in order to take corre- 
lations partially into account. The method usually im- 
proves the electronic-structure description, but it suffers 
from shortcomings associated with the [/-dependence of 
the calculated properties Unfortunately, there is usu- 
ally no obvious choice of the U value to be adopted; 
common choices are usually based either on experimen- 
tal input or are derived from constrained DFT calcula- 
tions, but neither of these approaches is entirely satis- 



factory. When dealing with phonon calculations, the U- 
dependence becomes even more critical since the phonon 
frequencies depend strongly on the unit-cell volume used 
for their evaluation, and the theoretical volume, in turn, 
depends on U. It is worth mentioning that even if an 
appropriate choice of U can accurately reproduce the 
binding energy of localized d states of transition-metal 
oxides, it is by no means guaranteed that the same U 
can accurately reproduce other properties of the same 
compound, such as the volumei^i2^ While most papers 
address the spin-phonon coupling by applying DFT-I-C/ 
methods, only a few deal with the dependence upon the 
U parameter 1^ In this paper, we will show that the spin- 
phonon coupling can strongly depend on the U parame- 
ter, and that such a dependence may give rise to artifi- 
cially large couplings. It is important to stress this mes- 
sage, often overlooked in the literature, in view of the 
increasing interest in ab-initio predictions of ferroelectric 
materials driven by spin-phonon coupling. 

In the last few years, another paradigmatic approach 
has been widely applied in solid-state materials science, 
namely the use of "hybrid functionals" that incorporate a 
weighted mixture of the exchange defined in the Hartree- 
Fock theory (but using the Kohn-Sham orbitals) and the 
density-functional exchange. The correlation term is re- 
tained from the density-functional framework. It is now 
widely accepted that the hybrid functionals outperform 
semilocal functionals, especially for bulk materials with 
band gaps»^"— It has also been shown that for low- 
dimensional systems such as semiconductor/oxides in- 
terfaces, the performance of hybrid functionals remains 
quite satisfactory.^*^ However, some doubts have very re- 
cently been put forward about the performance of hybrid 
functionals for certain structural configurations, e.g., sur- 
faces or nanostructures^^ 

While there is a plethora of different hybrid function- 
als, many of them defined empirically, a suitable func- 
tional derived on theoretical grounds is the so called 
FBEOj-SSiSa where the exchange mixing parameter as been 
fixed to one quarter as justified by a perturbation-theory 
calculation. A closely related functional, the Heyd- 
Scuseria-Ernzerhof (HSE) hybrid functional,— intro- 
duces yet another parameter n which splits the Coulomb 
interaction kernel into short- and long-range pieces while 
retaining the mixing only on the short-range component. 
It has been shown that this new hybrid functional;^ 
while preserving most of the improved performance of 
PBEO with respect to standard local and semi-local ex- 
change correlation functionals, greatly reduces the com- 
putational cost. For this reason, it is especially suitable 
for periodically extended systems, and is currently be- 
ing applied in many solid-state applications ranging from 
simple semiconductor systems to transition metals, lan- 
thanides, actinides, molecules at surfaces, diluted mag- 
netic semiconductors, and carbon nanostructures (for a 
recent review see Ref . Issl) . The HSE functional has been 
also used for phonon calculations for simple semiconduct- 
ing systems^iS or perovskite structures, ^'^'^^ where it 
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was shown the that phonon modes are much more accu- 
rately reproduced using the hybrid functionals than using 
GGA or LDA.^ 

Very recently, extended benchmarkings of the HSE 
method as well as other self-interaction-corrected ap- 
proaches have been presented for prototypical transi- 
tion metal oxides such as MnO, NiO and LaMnOa, 
and it has been shown that "HSE shows a remarkable 
quantitative agreement with experiments on most ex- 
amined properties"— and "HSE shows predictive power 
in describing exchange interactions in transition metal 
oxides." These recent studies further motivate us to 
use the HSE functional for our studies of the spin-phonon 
coupling; indeed, as we will see in Sec. IHI Bl the HSE- 
calculated phonon modes agree well with available exper- 
imental results for SrMnOa. The spin-phonon coupling 
effect, to the best of our knowledge, is totally unexplored 
by hybrid functional approaches. In this work we aim at 
filling this gap. 

Even when using the more efficient HSE functional, 
however, the use of hybrids entails an increased compu- 
tational cost which makes the calculation of phonon prop- 
erties of complex magnetic oxides very difficult. In this 
paper, we propose to circumvent this bottleneck by com- 
bining the HSE and DFT-fC/ approaches, i.e., choosing 
the appropriate U for each material by fitting to the HSE 
results for some appropriate materials properties. This 
provides a fairly efficient and affordable strategy that 
preserves the "HSE accuracy" for lattice constants, spin- 
phonon couplings, and related properties, while taking 
advantage of the computationally inexpensive DFT+U 
method for the detailed calculations. Further details will 
be given in Section HIl 

As far as the materials are concerned, it has been sug- 
gested that ^MnOs perovskites with A = Ca, Sr and 
Ba may be good candidates for spin-phonon-coupling 
driven multiferroicity.— i^i^ Furthermore, it has been re- 
ported that both SrMnOs and CaMnOa have a large 
spin-phonon coupling. ■^"■^^ We will show below that the 
spin-phonon coupling can depend strongly on the chosen 
U. We will also consider the simple perovskite LaMOa 
materials with M = Cr and Fe, which have Neel temper- 
atures above room temperature^ and large band gaps. 
We want to explore the possibility of using these two 
materials as building blocks for room temperature mul- 
tiferroics, e.g., in the form of double perovskites such as 
La2CrFe06. We have chosen these two classes of ma- 
terials in part for the reasons outlined above, but also 
because they are sufficiently "easy-to-calculate" for the 
benchmark and testing purposes of the present work, es- 
pecially when considering hybrid functionals. 

The paper is organized as follows. In Sec. [H] we re- 
port the computational details and describe our strategy 
for fitting U of the GGA+U calculations via a prelim- 
inary "exploratory" study using the HSE method. In 
Sec. mil we discuss: i) the effect of U on the frequency 
shift, in Sec. IHI Al n) the spin-phonon coupling effects in 
AMnOs with A = Ca, Sr, and Ba, in Sec. ilHBl in) the 



spin-phonon coupling effects in LaAf O3 with M = Cr, Fe, 
and the corresponding double perovskite La(Cr,Fe)03 in 
Sec. IHI CI and iv) prospects for the design of new multi- 
ferroics in Sec. IHI Dl In Sec. IIVI we give a summary and 
conclusions. Finally, in the Appendix, we give details 
about the methodology used here for a full HSE phonon 
calculation. 



II. METHODS 




FIG. 1: (Color online) ABO3 perovskite structure doubled 
along the [111] direction. A atoms (largest) shown in green; 
B atoms shown in violet; O atoms (smallest) shown in red. 

All of our calculations were performed in the frame- 
work of density-functional theory as implemented in 
the Vienna ab-initio simulation package (VASP-5.2)^i^ 
with a plane- wave cutoff of 500 eV. The DFT+U and 
HSE calculations were carried out using the same set of 
projector-augmented-wave (PAW) potentials to describe 
the electron-ion interactioni^SiSS The unit cell for sim- 
ulating the G-AFM magnetic ordering, where all first- 
neighbor spins are antiferromagnetically aligned, is dou- 
bled along the [111] direction as shown in Fig. [1] The 
same simulation cell is retained for the FM configura- 
tion in order to avoid numerical artefacts that could arise 
in comparing calculations with different effective k-point 
samplings. A 6 x 6 x 6 F-centered /c-point mesh is used. 

The phonon frequencies were calculated using the 
frozen-phonon method. Except for the case of 
La2CrFe06, the charge density and dynamical matrix re- 
tain the primitive 5-atom-cell periodicity for both the G- 
AFM and FM spin structures, so it is appropriate to an- 
alyze the phonons using symmetry labels from the prim- 
itive PmSm perovskite cell. The zone-center phonons 
decompose as ® F^ (plus acoustic modes), with the 
modes being polar. The zone-boundary modes at the 
R point (which appear at the F point in our 10-atom-cell 
calculations) decompose as © R^ © 2 R4 © R^ © R;^ , 
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FIG. 2; (Color online) Three idealized polar modes and _Rg mode in simple perovskites. (a) Slater mode; (b) Last mode; (c) 
Axe mode; (d) mode. 



where the three-fold degenerate is of special inter- 
est because it corresponds to rigid rotations and tilts of 
the oxygen octahedra. After each frozen-phonon calcu- 
lation of the zone-center phonons of our 10-atom cell, we 
analyze the modes to assign them according to these k- 
point and symmetry labels. For La2CrFe06, some of the 
modes mix (e.g., with R^); in these cases we report 
the dominant mode character. 

As for the polar phonons, these are often char- 
acterized in ABO3 perovskites in terms of three kinds 
of idealized modes as illustrated in Figs. [2^-c. The 
Slater mode (S-F4 ) of Fig. ^a) describes the vibration 
of B cations against the oxygen octahedra; the Last (L- 
F4 ) mode of Fig. ^b) expresses a vibration of the A 
cations against the BOq octahedra; and the Axe mode 
(A-r4 ) of Fig. [5Jc) represents the distortion of the oxy- 
gen octahedral^ The actual mode eigenvectors never be- 
have precisely like these idealized cases, but we find that 
they can be identified in practice as being mainly of 
one character, which is the one we report. These po- 
lar modes contribute to the low-frequency dielectric con- 
stant, and their softening in the high-symmetry paraelec- 
tric phase indicates the presence of a ferroelectric insta- 
bility. For the insulating compounds, we further calcu- 
lated their dielectric constants and Born effective charges 
using density-functional perturbation theory within the 
DFT+U context as implemented in VASP. The antifer- 
rodistortive (AFD) mode (i?^ mode), which describes 
the rotation oxygen octahedra, is also shown in Fig. [IJd) . 

We make use of the DFT+U method^ in the standard 
Dudarev implementation where the on-site Coulomb in- 
teraction for the localized 3d orbitals is parametrized by 
UoS = U — J (which we denote henceforth as simply 
using the PBEsol functional^ which has been shown to 
give a satisfactory description of solid-state equilibrium 
properties. We shall refer to this as PBEsol-l-(7. The lack 
of experimentally available data for our systems prevents 
us from extracting U directly from experiments; we will 
return to this delicate point shortly. 

The other functional we have used is HSEOG,^'' a 
screened hybrid functional introduced by Heyd, Scuse- 
ria, and Ernzerhof (HSE), where one quarter of the 
PBE short-range exchange is replaced by exact exchange, 
while the full PBE correlation energy is included. The 



range-separation parameter /i is set to fj, — 0.207A^^. 
The splitting of the Coulomb interaction into short- and 
long-range pieces, as done in HSE, allows for a faster 
numerical convergence with fc-points when dealing with 
solid-state systems. However, as previously mentioned, 
the application of the HSE approach to phonon calcula- 
tions for magnetic oxide systems remains very challenging 
in terms of computational workload. 

Scheme to fit U from hybrid calculations 

Here, we propose a practical scheme to perform rel- 
atively inexpensive DFT+U simulations that retain an 
accuracy comparable to HSE for the calculation of spin- 
phonon couplings. 

Our goal is to obtain a DFT+U approach that re- 
produces the dependence of the materials properties on 
the spin arrangement as obtained from HSE calculations. 
Naturally, the first and most basic property we would 
like to capture correctly is the energy itself. As it turns 
out, the energy is also a very sensible property on which 
one can base a [/-fitting scheme: The energy differences 
between spin configurations are directly related to the 
magnetic interactions or exchange constants, and these 
are known to depend on the on-site Coulomb repulsion 
U affecting the electrons of the magnetic species. (Typi- 
cally, in our compounds of interest, the value of U used in 
the simulations will play an important role in determin- 
ing the character of the top valence states; in turn, this 
will have a direct impact on the nature and magnitude 
of the exchange couplings between spins.) This depen- 
dence of the exchange constants on U makes this crite- 
rion a very convenient one for our purposes. Of course, 
such a fitting procedure does not guarantee our DFT+U 
scheme will reproduce correctly the phonon frequencies 
and frequency shifts between different spin arrangements 
obtained from HSE calculations. In that sense, we are 
relying on the physical soundness of the Hubbard- C/ cor- 
rection to DFT; as we will see below, the results are quite 
convincing. 

Obtaining U from the energy differences has an addi- 
tional advantage: It allows us to devise a very simple fit- 
ting procedure that relies on a minimal number of HSE 
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FIG. 3: (Color online) Variation of PBEsol+t/ total energy 
difference AE = Eafm - Efm with U for LaCrOs, Aj5+^ 
is the energy difference with volume fixed to the optimized 
volume from HSE and A_B^^ relaxed volume for each U. (The 
fitted value of U occurs at the crossing with the AE^^^ value.) 



calculations. In essence, these are the steps we follow. 
We relax the cubic structure separately for G-AFM and 
FM spin configurations using HSE and obtain the total- 
energy difference A£:(HSE) =£:afm (HSE) - £:fm(HSE). 
We then carry out a series of PBEsol+ U calculations in 
which U is varied from to 6eV and obtain AE{+U) ~ 
EAFui+U) — Epm{+U) for each U. The U is then chose 
such that AE{+U)^ A£;(HSE). 

Figure [3] shows the results of our fitting for LaCrOs 
by using two slightly different approaches. In one case 
we ran the PBEsol-|- U simulations at the HSE-optimized 
volumes, and in the other case we performed volume re- 
laxations at the PBEsol-|- U level for each U value consid- 
ered. It can be seen that such a choice does not have a big 
effect on the computed frequency shifts, and both result 
in the same U value. Hence, the U values reported here 
were obtained by running PBEsol-|- U simulations at the 
HSE volumesii^ We obtained 3.0, 2.8, 2.7, 3.8, and 
5.1 eV, respectively, for CaMnOa, SrMnOs, BaMnOa, 
LaCrOa and LaFeOai^i We then do detailed phonon cal- 
culations using this value of U in the PBEsol-t-L^ cal- 
culations to investigate spin-phonon coupling effects in 
AMnOa and LaMOs. 

As we will show in Sec. IIIIl we have tested this pro- 
posed scheme and found that it works quite well. In par- 
ticular, the phonon frequency shifts Aw = wafm ~ ^^fm 
computed using PBEsol-|- U with the fitted U are almost 
the same as those obtained using HSE directly. Since a 
previous study has concluded that HSE works "remark- 
ably well" for transition-metal oxideS)^ we believe this 
approach can be used with confidence. 

Finally, we note that the direct HSE calculations can 
be rather heavy, even though we only have 10-atom cells; 
the presence of magnetic order and the need to calculate 
the phonons and the spin-phonon couplings makes the 
calculations challenging.^^ We circumvent this difficulty 
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FIG. 4: (Color online) Effect of U on spin-phonon coupling 
in SrMnOs. (a) FJ and phonon frequency shifts (Aoj — 
ijJAFM — i^fm) versus U. (b) S-F^ and phonon frequencies 
versus U for AFM and FM states. 

by splitting a single run of HSE frozen-phonon calcula- 
tions into several parallel runs, each one computing the 
forces for symmetry-independent atomic displacements. 
We then use the forces calculated in these runs to con- 
struct the force-constant matrix, and by diagonalizing 
this matrix, we obtain the phonon frequencies and eigen- 
vectors. Further details of this procedure are presented 
in the Appendix. 



III. RESULTS AND DISCUSSION 
A. Effect of U on the frequency shifts 

Let us begin by showing some representative results 
of the U -dependence of the phonon frequencies and fre- 
quency shifts for the ^MnOa and LaA/Os compounds. 
Here we focus on the case of SrMnOa, which has been 
predicted to exhibit a large spin-phonon coupling based 
on GGA+U simulations using 1.7 eV; more precisely, 
a giant frequency shift Auj= 230 cm^^ has been reported 
for the Slater mode at the theoretical equilibrium state»^ 

We computed the phonon frequencies and frequency 
shifts of the ideal cubic perovskite phase of SrMnOa us- 
ing U values in the range between and 8eV. Figure 0] 
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TABLE I: Lattice constant a, local Mn magnetic moment m, 
and band gap -Egap (zero if blank) for CaMnOa, SrMnOa and 
BaMnOa from PBEsol+f/ (+17) and HSE methods. 







a (A) 


m (^b) 




CaMnOg-AFM 


+U 
HSE 


3.73 
3.73 


2.79 
2.81 


0.35 
1.70 


CaMnOs-FM 


+U 
HSE 


3.74 
3.73 


2.84 
2.81 




SrMnOa-AFM 


+U 
HSE 


3.80 
3.80 


2.80 
2.85 


0.40 
1.45 


SrMnOa-FM 


+U 
HSE 


3.81 
3.81 


2.89 
2.91 


0.00" 


BaMnOa-AFM 


+U 
HSE 


3.91 
3.91 


2.86 
2.92 


0.13 
1.20 


BaMnOa-FM 


+U 
HSE 


3.93 
3.93 


3.14 
3.13 





"Half-metallic. 



TABLE II: Phonon frequencies and frequency shifts calcu- 
lated for CaMnOa, in cm~^. Frequencies calculated using 
HSE and FBEsol-t-?/: uj^^^ and Frequency difference 

between HSE and PBEsol-|-[/: = -uj^^^. Frequency 
shift: Aoj = UAYM — ojym- Difference of frequency shifts be- 
tween HSE and PBEsol+17: Ao;^ = Aa;+^ - Aa;^^'^. 





HSE 
'^AFM 


^AFM 


HSE 
^FM 


^FM 


A HSE 




L-F4 


126 


-29 


-81 


-18 


207 


-10 


S-F4 


272 


-11 


186 


-3 


86 


-8 




241 


-48 


183 


-72 


58 


24 




658 


10 


636 


11 


22 


-1 


R4 


430 


-46 


413 


-49 


17 


2 




-204 


-18 


-218 


-20 


14 


2 


A-F4 


569 


35 


557 


40 


12 


-5 




437 


-11 


448 


-14 


-11 


3 


R2 


890 


-27 


885 


-31 


5 


4 


R4- 


160 


-20 


157 


-21 


3 


1 



shows our results for the key modes that determine the 
occurrence of ferroelectricity, i.e., the F4 polar phonons 
and the R^ antiferrodistortive (AFD) mode. The Rg" 
mode involves antiphase rotations of the Og octahedra 
around the three principal axes of the perovskite lattice; 
such a mode is soft in many cubic perovskites, and it of- 
ten competes with the FE instabilities to determine the 
nature of the low-symmetry phases. 

From Fig. |4^ we see that the r4 and Rg" modes de- 
pend strongly on the chosen U . Notably, the frequency 
shift Aw of the S-r4 mode can change from 400 cm~^ 
to cm^^ as U increases from to 8eV. The Aw of the 
A-r4 and R^ modes also depend on the t/, while the 
L-F^ mode is nearly insensitive to U within this range 
of values. Interestingly, the large frequency shift of S-r4 
is related to the strong {/-dependence of the FM S-F4 
mode, as shown in Fig. ^p. 

Hence, our results show that the magnitude of the 
spin-phonon coupling has a significant dependence on the 
value of U employed in a DFT-|-[/ calculation, and that 
such a dependence is particularly strong for some of the 
key modes determining the structural (FE/AFD) insta- 
bilities of cubic perovskite oxides. It is thus clear that 
choosing an appropriate U is of critical importance if we 
want to avoid artificially "strong" couplings. 



B. Spin-phonon coupling in CaMnOs, SrMnOa, 
and BaMnO:^ 

We studied the cubic phase of CaMnOs, SrMnOa, 
and BaMnOs using the PBEsol+J7 approximation with 
U values for Mn's 3d electrons (i.e., U = 3.0 eV for 
CaMnOa, U = 2.8 eV for SrMnOa, and U = 2.7 eV 
for BaMnOa) that were determined as described in Sec- 
tion |TT1 The basic properties that we obtained are listed 
in Table HI Our PBEsol-t-fJ method gives good lattice 



constants and local Mn magnetic moments compared 
with HSE. As regards the metallic or insulating char- 
acter, the PBEsol-|-(7 agrees qualitatively with the HSE 
result in most cases, although the band gap is under- 
estimated. In FM-SrMnOa there is a clear discrepancy 
between PBEsol+f/ and HSE: the HSE calculations pre- 
dict a half-metallic state, while we obtain a metal with 
PBEsol-K[7. 

We have also calculated the Born effective charges 
for insulating AFM configurations of CaMnOa (SrMnOa, 
BaMnOa). We obtain Zca/Sr/Ba=2.60 e (2.58 e, 2.73 e), 
ZM„=7.35e (7.83 e, 9.42 e), Zo||=-6.55e (-6.93 e, 
-8.07e), Zo^=-1.70e (-1.74e, -2.04e), where e is the 
elementary charge and and Zq^ refer, respectively, 
to the dynamical charges defined for an atomic displace- 
ment parallel and perpendicular to the Mn-0 bond. The 
anomalously large Zmu and Zq^^ charges of AFM state 
of CaMnOa, SrMnOa and BaMnOa are strongly reminis- 
cent of the results obtained for ferroelectric perovskite 
oxidesr^ and suggest the possible proximity of a polar 
instability that might be triggered by an appropriate ex- 
ternal (e.g., strain) fieldj^i^SiS^ 

The phonons at F and R for different spin orders for 
CaMnOa, SrMnOa and BaMnOa are shown in Tables [Hi 
mil and IIVI respectively; the phonon modes are ordered 
by descending Aw^^^. The first thing to note from the 
table is that the frequencies obtained for SrMnOa 's po- 
lar modes (i.e., 177, 187 and 494 cm~^) agree well with 
available low temperature experiment results (i.e., 162, 
188 and 498 cm^^ taken from Ref.[73). Thus, our results 
provide additional evidence that the HSE scheme renders 
accurate phonon frequencies for magnetic oxides. The 
second thing to note from these tables is that the phonon- 
frequency shifts (Aw) obtained with PBEsol-|-I7 are in 
overall good agreement with the HSE results, which sug- 
gests that our strategy to fit the value of C/ is a good 
one. 
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TABLE III: Calculated phonon frequencies and frequency TABLE IV: Calculated phonon frequencies and frequency 
shifts, as in Table |II1 but for SrMnOa. shifts, as in Table |IT1 but for BaMnOa. 





, ,HSE 
'^AFM 


^AFM 


, ,HSE 
"^FM 


^FM 


^^HSE 






, ,HSE 
^AFM 


/ A 

^AFM 


, ,HSE 
■^FM 


/ A 

^FM 




Alu^ 


s-r4 


177 


40 


-96 


-16 


273 


56 


s-r4 


-274 


117 


-369 


227 


95 


-110 




64 


-135 


-6 


-121 


70 


-14 


R2 


670 


-28 


631 


-59 


39 


31 




292 


-38 


253 


-54 


39 


16 


^5 


317 


-34 


281 


-28 


36 


-6 




81 1 


zu 


7SQ 
( oy 




99 


_1 


R + 


ouo 


—1 
10 


000 


-99 
zz 


98 


Q 




571 


12 


553 


2 


18 


10 


R5: 


221 


-28 


200 


-67 


21 


39 


R4 


424 


-47 


410 


-62 


14 


15 


R4 


403 


-48 


384 


-60 


19 


12 


L-r4 


187 


-16 


177 


-7 


10 


-9 


A-r4 


423 


-12 


404 


11 


19 


-22 


A-r4 


494 


22 


487 


17 


7 


5 


^3 


423 


16 


407 


-20 


16 


36 




412 


-12 


409 


-12 


3 







200 


-6 


195 


-6 


5 





R4 


159 


-6 


157 


-7 


2 





R4 


150 


-5 


147 


-6 


3 


1 



Third, our results offer information about the depen- 
dence of the structural instabilities on the choice of the 
A-site cation and on the magnetic arrangement. We find 
that the AFD instability {R^ mode) becomes weaker and 
disappears, for both AFM and FM states, as the size of 
the A-site cation increases (The effective ionic radii esti- 
mated by Shannon^^ are rca = 1.34 A, rsr = 1.44 A, and 
TBa = 1.61 A.) In contrast, the ferroelectric instability 
(Slater mode) becomes stronger as the A-site cation be- 
comes bigger. This is the usual behavior that one would 
expect for these two instabilities of the cubic perovskite 
structure, and has been recently examined in detail by 
other authors.— 1^^'^° It clearly suggests that some of these 
compounds could display a magnetically ordered ferro- 
electric ground state. In particular, this could be the 
case for SrMnOa and BaMnOsi in these compounds, the 
AFD instability becomes weaker or even disappears, so 
that it can no longer compete with and suppress the FE 
soft mode. 

The phonon frequency shift Alu between the AFM 
and FM magnetic orders is shown in Fig. \5\ We 



find that the Slater modes exhibit a very consider- 
able spin-phonon coupling for all three compounds, with 
Aw > 100 cm^^. A sizable effect is also obtained for 
the mode (Acj ^ 50 cm~^) in all cases. Additionally, 
the Last mode shows a very large effect for CaMnOa, 
and the spin-phonon coupling for R^ is significant in the 
case of SrMnOs. The coupling is relatively small, or even 
negligible, for all other modes. 

We will not attempt any detailed interpretation of our 
quantitative results here, but some words are in order. 
In the recent literature, it has often been claimed that 



(a) 




1 1 1 1 r— : — t: — ~i r 



CaMnOg — Q- 
SrMnOg — ©- 
BaMnOg — 0- 




^"^4 ^5 L-r4 A-r4 R3 R4 R4 Rg"^ R2 R5 



FIG. 5: (Color online) Frequency shifts Aoj — ojafm — (jJfm 
in ^MOa from HSE calculations. 



(b) 



FIG. 6: (Color online) Changes to metal-O-metal bond angle 
(solid blue line) resulting from Slater (a) and (b) phonon 
modes. Open and solid dots indicate ideal and displaced po- 
sitions respectively. Small red dots are oxygen; larger purple 
dots are metal atoms. 
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TABLE V: Lattice constant a, local magnetic moment m, band gap -Egap, dielectric constant e, and Born effective charge Z 
for LaCrOa, LaFeOs and LaaCrFeOe (LCFO) from GGA+U (+[/) and HSE methods. 







a (A) 


m (un) 


■E^ean (eV) 


£00 


£0 


Zj,p, (e) 


Znr/Vr (6) 


Zon (e) 


Zo^ (e) 


LaCrOa-AFM 


+U 
HSE 


3.87 
3.87 


2.86 
2.75 


2.01 
3.06 


6.12 


93.71 


4.53 


3.34 


-3.50 


-2.19 


LaCrOa-FM 


+U 
HSE 


3.88 
3.88 


2.84 
2.74 


1.70 
2.49 


6.08 


236.19 


4.52 


3.61 


-3.62 


-2.26 


LaFeOs-AFM 


+U 
HSE 


3.90 
3.91 


4.19 
4.10 


2.20 
3.25 


6.37 




4.47 


3.75 


-3.41 


-2.43 


LaFeOg-FM 


+U 


3.91 


4.38 


1.57 


6.19 




4.48 


4.02 


-3.59 


-2.46 


HSE 


3.93 


4.26 


2.20 














LCFO-AFM 


+U 
HSE 


3.88 
3.88 


2.71/4.30 
2.63/4.19 


1.98 
2.78 


6.35 


713.57 


4.51 


3.02/4.22 


-3.50 


-2.33 


LCFO-FM 


+U 
HSE 


3.89 
3.89 


3.01/4.28 
2.89/4.18 


2.01 
3.04 


6.18 


262.40 


4.49 


3.21/3.87 


-3.45 


-2.30 



AFD distortions are expected to couple strongly with 
the magnetic structure of perovskite oxides, as they con- 
trol the metal-oxygen-metal angles that are critical to 
determine the magnitude of the magnetic interactions 
that are dominant in insulators.— [The hopping param- 
eters between oxygen (2p) and transition metal (3c?) or- 
bitals are strongly dependent on such angles, as are the 
effective hoppings between 3c? orbitals of neighboring 
transition metals. As we know form the Goodenough- 
Kanamori- Anderson rules, the super-exchange interac- 
tions are strongly dependent on such hoppings.] In prin- 
ciple, such an effect should be reflected in our computed 
frequency shifts; for such R^^ modes we obtain Aw values 
ranging between 10 cm~^ and 70 cm~^. Interestingly, fol- 
lowing the same argument, one would conclude that the 
Slater and F^ modes should have a similar impact on 
the magnetic couplings in these materials, as both involve 
changes in the metal-oxygen-metal angles due to the rel- 
ative displacement of metal and oxygen atoms, as shown 
in Fig.ini Further, the Slater modes also affect the metal- 
oxygen distances, as they involve a significant shortening 
of some metal-oxygen bonds (Fig.|6^). Hence, from this 
perspective, the large Acj values obtained in our calcula- 
tions for these two types of modes are hardly surprising. 
Finally, the result obtained for the Last mode in the case 
of CaMnOs (i.e., Aui w 200 cm"-"^) is clearly anomalous 
and unexpected according to the above arguments. We 
speculate that it may be related to a significant Ca-0 
interaction that alters the usual super-exchange mecha- 
nism and favors a FM interaction. While unusual, effects 
that seem similar have been reported previously for other 
compounds 

C. Spin-phonon coupling in LaCrOs, LaFeOs, and 
La2CrFe06 

In the previous Section we have shown how spin- 
phonon coupling effects can crucially depend on the na- 
ture of the A-site cations in AMnOa compounds. Now 



we will focus on the change of B-site cation to investigate 
the spin-phonon couplings in the LaMOa compounds 
(M=Cr, Fe) and the double perovskite La2CrFe06. For 
our PBEsol-|-/7 calculations we used U = 3.8 eV for Cr 
and U = 5.1 eV for Fe, which were determined as de- 
scribed in Section ini We obtained these C/'s from calcu- 
lations for LaCrOa and LaFeOa and used the same values 
for the PBEsol-hf/ study of La2CrFe06. 

The basic computed properties of LaCrOa, LaFeOa 
and La2CrFe06 are presented in Table fVl In all cases, we 
obtain insulating solutions for AFM and FM spin orders, 
both from HSE and PBEsol-t-[/ calculations. Therefore, 
we also calculated the Born effective charges and optical 
dielectric constants within PBEsol+J7 by using density- 
functional perturbation theory as implemented in VASP. 
The static dielectric constants were also calculated for 
compounds which do not have an unstable polar mode. 
It is evident that and Zq^ are essentially identical 
for the three compounds and insensitive to the spin or- 
der. On the other hand, the Zci/fc a-nd charges 
of LaCrOa and LaFeOa increase in magnitude when the 
spin arrangement changes from AFM to FM, and de- 



TABLE VI: Calculated phonon frequencies and frequency 
shifts, as in Table HIl but for LaCrOs. 





, ,HSE 
'^AFM 




, ,HSE 


"^FM 


A^HSE 






230 


-24 


178 


-35 


52 


11 




361 


-26 


317 


-27 


44 


1 


L-r4 


69 





48 


-3 


21 


3 


R4 


389 


-29 


372 


-30 


17 


1 




414 


2 


430 





-16 


2 


A-F4 


659 


32 


644 


30 


15 


2 




-213 


1 


-228 


1 


15 







683 


5 


669 


5 


14 





R4 


89 


-4 


81 


-4 


8 





R2" 


840 


-17 


843 


-19 


-3 


2 
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crease in the case of the double perovskite La2CrFe06. 
Table |V] also shows that the optical dielectric constants 
£00 are very close to 6 for these three materials, and are 
independent of the magnetic order. However, the static 
dielectric constants Eq of LaCrOa and La2CrFe06 change 
very significantly when moving from AFM to FM. Also, 
the static dielectric constant of AFM-La2CrFe06 is very 
large due to the very small frequency of the Last phonon 
mode. The static dielectric constants of LaFeOa are not 
shown in Table IVl because the Last modes are unstable, 
and strictly speaking they are not well defined. (Roughly 
speaking, in all such cases we would have £0 — ^ 00, as the 
cubic phase is unstable with respect to a polar distor- 
tion.) 

In Tables IVII and I VIII we show the phonon frequen- 
cies as calculated at the HSE and PBEsol-l- U levels for 
LaCrOs and LaFeOa respectively. Further, Table IVIIII 
shows the PBEsol-|- (7 results for La2CrFe06. As was the 
case for the AMnOa compounds of the previous section, 
we find here as well that the AFM-FM frequency shifts 
computed with our PBEsol+ (7 scheme reproduce well the 
HSE resuhs. 

Our results also show that the phonons of the LaMOs 
compounds exhibit some features that differ from those 
of the AMnOa materials. First, the octahedral rotation 
mode (R^) is unstable for all the La-based compounds, 
and it is largely insensitive to the nature of the B-site 
cation. Second, the energetics of the FE modes is very 
different. In the case of the Mn-based compounds, the 
Slater mode is the lowest-frequency mode, and in some 
cases it becomes unstable, thus inducing a polarization, 
just as the unstable Slater mode induces ferroelectricity 
in BaTiOa.^*' However, in the La compounds the lowest- 
frequency mode is the Last phonon mode, and in the 
cases of LaFeOa and La2CrFe06 this Last mode is un- 
stable and might induce ferroelectricity; this situation is 
more analogous to what occurs in FbTiOg^ or BiFeOa. 

The spin-phonon coupling effects computed for the 
LaMOa compounds are given in Tables lVTllVIH and lVIIl] 
and are summarized in Fig. [71 Here, the first thing to 
note is that the magnitude of the effects is significantly 
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FIG. 7: (Color online) Frequency shifts for LaMOa (HSE for 
LaCrOs and LaFeOs; PBEsol+[/ for La2CrFe06). 



smaller than for the LaA/Oa compounds (note the dif- 
ferent scales of Figs. [5] and [7]). Second, we observe that 
the largest effects are associated with the Fj^ and Slater 
modes (with Aw « 45 cm~^ for LaCrOa). This is consis- 
tent with the point made above that these modes disrupt 
the metal-oxygen-metal super-exchange paths. In com- 
parison, in this case we obtain a relatively small effect for 
the R^ modes, which show Aw values (< 20 cm~^) that 
are comparable to those computed for most of the phonon 
modes considered. Finally, for LaCrOa and LaFeOa we 
observe that the F phonon frequencies decrease as the 
spin order changes from AFM to FM, in line with what 
was observed for the AMnOa compounds. In contrast, 
the r modes of double-perovskite La2CrFe06 increase in 
frequency when the spin order changes to FM. 



D. Prospects for the design of new multiferroics 

Let us now discuss several possible implications of our 
results regarding the design of novel multiferroic materi- 



TABLE VII: Calculated phonon frequencies and frequency TABLE VIII: Calculated phonon frequencies and frequency 
shifts, as in TableHl but for LaFeOs. shifts, as in TableHIl but for La2CrFe06 from PBEsol+C/.^ 





, ,HSE 
^AFM 


"^AFM 


Wpjvi 




AwHSE 






(^AFM 


OJpM 


Aw 


R2 


796 


-38 
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1. BaMnOs-based materials 

As already mentioned in the Introduction, the spin- 
phonon couphng can trigger multiferroicity in some ma- 
terials if, by application of an epitaxial strain or other 
perturbation, a FM-FE state can be stabilized with re- 
spect to the AFM-PE ground state. Indeed, the three 
investigated AMnOs compounds have been suggested as 
candidates to exhibit strain-induced multiferroicity by 
three research groups. ^■•^'^■^'^ According to our calcula- 
tions, we can tentatively suggest that BaMnOs might 
show multiferroicity even without strain applied, pro- 
vided the material can be grown as a distorted perovskite 
(the most stable polymorph of BaMnOa adopts instead a 
structure with face-sharing octahedra^^'^**). Note that it 
may be possible to enhance the stability of BaMnOa's cu- 
bic perovskite phase by partial substitution of Ba by Ca 
or Sr. In fact, ideally one would try to obtain samples of 
Bai_a;Sra;Mn03 or Bai-^^CaajMnOa with x large enough 
to stabilize the perovskite phase, and small enough for 
the FE instability associated to the FM Slater mode to 
dominate over the R^" and AFM Slater modes. Appar- 
ently this strategy to obtain new multiferroics has re- 
cently been realized experimentally by Sakai et alr^ in 
Bai_a;Sra;Mn03 solid solutions. An alternative would be 
to consider CaMnOs/BaMnOs or SrMnOs/BaMnOa su- 
per lattices with varying ratios of the pure compounds, 
and perhaps tuning the misfit strain via the choice of 
substrate. 

Another intriguing possibility pertains to the mag- 
netic response of the AMnOs compounds that display 
an AFM-PE ground state and a dominant polar insta- 
bility of their FM phase. Again, this could be the case 
for some Bai_2;Sr2;Mn03 and Bai-a^CaajMnOs solid so- 
lutions with an appropriate choice of x. By applying a 
magnetic field to such compounds, it might be possible 
to switch them from the AFM ground state to a FM 
spin configuration and, as a result, induce a ferroelectric 
polarization. 



2. Double-perovskite La2 CrFeOe 

The double perovskite La2CrFe06 has been intensively 
studied to examine its possible magnetic order through 
— 3d^ superexchange. However, its magnetic ground 
state has long been debated. Pickett et al.^^ predicted 
that the ferrimagnetic (FiM) ground state with a net 
spin moment of 2/iB/f.u. is more stable than the FM 
one with T/ie/f.u. (This FiM order can be viewed 
as a G-AFM configuration in which, for example, all 
Fe spins are pointing up and all Cr spins are point- 
ing down.) However, from GGA and LDA-|-[/ calcu- 
lations, Miura et al. found that the ground-state mag- 
netic ordering of La2CrFe06 is FiM in GGA, but that 



even a small U in LDA+U makes it FM42 The exper- 
imental picture is also unclear. Ueda et al. have grown 
a (lll)-oriented LaCrOs/LaFeOs superlattice which ex- 
hibits FM ordering, although the measured saturation 
magnetization is much smaller than expected.'*^ Very re- 
cently, Chakraverty et al. reported epitaxial La2CrFe06 
double-perovskite films grown by pulsed-laser deposition, 
and their sample exhibits FiM with a saturation magne- 
tization of 2.0 ± 0.15 AiB/f.u. at 



Our HSE calculations for La2CrFe06 with the atomic 
(oxygen) positions relaxed in the cubic structure shows 
the magnetic ground state has a FiM spin pattern lead- 
ing to a FiM structure with a net magnetization of 
l.Se^iB/f.u. This is consistent with LDA^^ and GGA§1 
calculations, as well as being heuristically consistent with 
the experimental report of AFM ordering in La2CrFe06 
solid-solution films-SS, However, our HSE calculation 
shows that the energy difference between the FiM and 
FM states is very small: FiM is only 0.8meV/f.u. lower 
in energy than FM. In addition, the GGA+U with U fit- 
ted to LaCrOs and LaFeOa {U —3.8 and 5.1 eV for Cr and 
Fe, respectively) results in a FM magnetic ground state 
having a total energy 34.9 meV lower than that of the 
FiM. By doing a fitting of the U parameters for Cr and 
Fe directly to i?AFM — £^fm of La2CrFe06 as computed by 
PBEsol-hC/ and HSE, instead of for LaCrOa and LaFeOg 
separately, we find parameters of C/=3.0 and 4.1 eV for 
Cr and Fe respectively (the phonon properties predicted 
from these are very close to our previous results). Using 
these parameters, we find that the FiM ground state is 
0.7 meV lower in energy than the FM state. For compar- 
ison, a straight PBEsol calculation (with U — 0) predicts 
that the energy of the FiM ground state is 596meV/f.u. 
lower than that of the FM state. Clearly, U should be 
chosen carefully in order to obtain the correct ground 
state of La2CrFe06. 



According to our HSE calculations, the magnetic 
ground state of La2CrFe06 is FiM, but with the FM 
state lying only very slightly higher in energy. This may 
be the reason why questions about the magnetic ground 
state of La2CrFe06 have long been debated; the energy 
diff'erence is so small that external perturbations (e.g., 
the epitaxial strain in a superlattice^) or variations in 
U between different LDA-I-C/®'' and GGA+U calculations 
may bring the FM energy below that of the AFM. Taken 
together with our results, shown in Table IVHll that the 
Last mode is close to going soft in this material, this sug- 
gests that La2CrFe06 might be a good candidate for a 
material in which multiferroic phase transitions could be 
induced, similar to what was shown for SrMnOa-^ and 
SrCoOs.^^ Because the spin ordering is so delicate, it 
seems likely that a small misfit strain could be enough to 
trigger such a transition. 
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IV. SUMMARY AND CONCLUSIONS 

In this work, we have studied the spin-phonon coupling 
for transition-metal oxides within density-functional the- 
ory. From the computational point of view, an accurate 
description of the electronic, structural, and vibrational 
properties on an equal footing is a prerequisite for a re- 
liable study of the coupling between spins and phonon. 
Taking note of the increasing evidence that hybrid func- 
tionals are suitable for this task, we have adopted the 
Heyd-Scuseria-Ernzerhof (HSE) screened hybrid func- 
tional for the present work. However, the accuracy of 
the HSE results comes at the cost of an increase of com- 
putational load, so that a full frozen-phonon calculation 
of the phonon modes remains prohibitively expensive in 
many cases. We propose to overcome this limitation by 
carrying out calculations at the DFT+U level using U 
parameters that have been fitted to HSE results for total- 
energy differences between spin configurations. Our re- 
sults show that the resulting DFT-I-J7 scheme reproduces 
the HSE results very accurately, especially in regard to 
the spin-phonon couplings of interest here. 

As regards the direct HSE phonon calculations, we 
have developed an approach in which we split the calcu- 
lation into separate, simultaneous frozen-phonon calcula- 
tions for different symmetry-adapted displacement pat- 
terns, and then combine the results to calculate and di- 
agonalize the dynamical matrix, thus accelerating these 
calculations significantly. 

Our important results can be summarized as follows. 
First, we have shown that the choice of C/ is a big concern 
in such studies, since the spin-phonon coupling can de- 
pend very strongly on U. As an alternative to extracting 
U from experimental studies, we propose here to obtain it 
by fitting to HSE calculations as illustrated above. Sec- 
ond, we have studied CaMnOa, SrMnOa and BaMnOa, 
focusing on trends in the spin-phonon coupling due to 
the increase of the A cation size. Based on the strain 
couplings and the spin-phonon interactions, we suggest 
theoretically that BaMnOa is more likely to show fer- 
roelectricity under tensile strain and, furthermore, that 
A-site substitution by a cation with smaller size may in- 
duce multiferroicity even without external strain. Third, 
we find that in the AMnOs materials class with A=Ca., 
Sr, and Ba, the frequency shift decreases as the A cation 
radius increases for the F phonons, while it increases for 
R-point phonons. Fourth, we have shown that chang- 
ing B-site cations may also have important effects on the 
dielectric properties: in LaMOa with M=Cr, Fe, and 
Cr/Fe, the phonon frequencies at F decrease as the spin 
order changes from AFM to FM for LaCrOa and LaFeOa, 
but they increase for the double perovskite La2(CrFe)06. 
Finally, we have shown that the polar phonon modes of 
the investigated perovskites tend to display the largest 
spin-phonon couplings, while modes involving rotations 
of the Og octahedra present considerable, but generally 
smaller, effects. Such observations may be relevant as 
regards current efforts to obtain large magnetostructural 



(and magnetoelectric) effects. 

We hope that our study will stimulate further work 
leading to rational design and strain engineering of mul- 
tiferroicity using spin-phonon couplings. 

Appendix: Efficient phonon calculations with hybrid 
functionals 

Even though we only have a 10-atom cell, we found 
that it can be quite expensive to use the HSE functional 
to carry out the needed spin-polarized calculations of 
phonon properties We overcome this limitation as fol- 
lows. First, we use symmetry to limit ourselves to sets 
of displacements that will block-diagonalize the force- 
constant matrix. For example, for the polar modes we 
move the cations along x and the O atoms along x and 
y. We then carry out self-consistent calculations on these 
displaced geometries, and from the forces we construct 
the relevant block of the force-constant matrix. Second, 
while the standard VASP implementation uses a "central 
difference" method in which ions are displaced by small 
amounts in both positive and negative displacements, we 
save some further effort by displacing only in the posi- 
tive direction. Finally, we note that the forces resulting 
from each pattern of atomic displacements can be calcu- 
lated independently, allowing us to split the calculation 
in parallel across independent groups of processors and 
thus further reduce the wall-clock time. 

We have checked the accuracy of this approach for FM- 
LaCrOa using PBEsol+C/ (C/=3.8eV) and for G-AFM- 
BaMnOa using HSE. For each case, we compared the re- 
sults of the standard implementation of the VASP frozen- 
ion calculation of phonon frequencies with the revised ap- 
proach described above. We take the ion displacements 
to be 0.015 A in all cases. We find that the RMS er- 
ror of ten different phonon frequencies is 1.2 cm^^ for 
PBEsol-l-C/ and 7.3 cm~^ for HSE. These results suggest 
that this method has acceptable accuracy with reduced 
computational cost. We propose that it could be used 
also for cases of lower symmetry and larger cells, thus 
making the HSE phonon calculations at F affordable in 
general. 

Recently, a revised Perdew-Burke-Ernzerhof func- 
tional base on PBEsol, which we refer to as HSEsol, was 
designed to yield more accurate equilibrium properties 
for solids and their surfaces. Compared to HSE, signifi- 
cant improvements were found for lattice constants and 
atomization energies of solidsi^ We also checked the ef- 
fect of using HSEsol on our calculations, as shown in 
Table HX] From this table, we can see that the phonons 
calculated with HSEsol are very close to those from HSE. 
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